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We rigorously establish the asymptotic equivalence between the height function of interacting 
dimers on the square lattice and the massless Gaussian free field. Our theorem explains the micro¬ 
scopic origin of the sine-Gordon field theory description away from the free fermion point, which has 
previously been elusive. We use a novel technique, based on the combination of discrete holomor- 
phicity with exact, constructive, renormalization group methods, which has the potential of being 
applicable to a variety of other non-integrable models at or close to criticality. 
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High temperature superconductivity and the physics 
of Resonance Valence Bonds (RVB) [1] was the original 
motivation for studying two-dimensional (2D) quantum 
dimers, which later became an important model for frus¬ 
trated magnetism, cold bosons, and many other systems 
with hard constraints [2]. In these contexts also classical 
dimers are of interest, not only because they capture the 
high temperature physics of their quantum counterpart, 
but also because for special values of the parameters the 
quantum static correlations can be expressed in terms 
of the classical ones [3]. The properties of a wide class 
of classical dimer models can be understood by using a 
celebrated result of half a century ago, the Kasteleyn the¬ 
orem [4], ensuring exact solvability and explicit expres¬ 
sions of the correlations, which can be written in terms of 
Pfaffians. By using this result and the above mentioned 
equivalence, the correlations of certain quantum dimer 
models at special values of the parameters on the square 
[3] and triangular lattice [5] were computed, finding a 
power law (critical), and an exponential (massive) large 
distance decay, respectively. 

However, exact solvability is limited to a special class 
of systems, and further progress in our understanding of 
the physics of dimers requires the analysis of what hap¬ 
pens away from integrability. We consider a prototypical 
non-solvable dimer model obtained by assuming a local 
interaction between parallel dimers: given a periodic box 
A C of side L (with L even), the partition function is 


za(a,to)= [ n 
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where A^a is the set of dimer coverings oi A, X = v/T 
with T the temperature, P is a plaquette (face of Z^) 
and Np{M) = 1 if the plaquette P is occupied by 
two parallel dimers in M, and Np{M) = 0 otherwise; 
the TO-dependence in the reference weight = 

1 -I- (5jqm(— l)®i tunes the distance from criticality; A 
tunes the distance from solvability, with A > 0 corre¬ 
sponding to a local attractive interaction. This model, 
in the m = 0 case, describes polar crystals [6] and it 
was recently reconsidered in [7-12] where its connection 


with quantum dimer models, RVB physics and large spin 
quantum anti-ferromagnets was worked out in detail and 
used to infer informations on the RVB spin-liquid order 
parameters. MonteCarlo simulations show the presence 
of non-universal anomalous exponents in the dimer cor¬ 
relations decay. This confirms the general picture that 
the asymptotic properties can be captured by a quantum 
field theory (QFT) of the sine-Gordon type, the funda¬ 
mental field being a coarse-grained version of the height 
function. Using this effective QFT description, several 
informations were derived about the phase diagram, in¬ 
cluding the Kosterlitz-Thouless universality of the phase 
transition from a liquid to a crystalline phase. The same 
effective description is believed to be applicable to a va¬ 
riety of dimer and interface models, and it is at the basis 
of our current understanding of their physics. However, 
while the validity of the QFT description is supported a 
posteriori by the agreement of its prediction with simu¬ 
lations, a purely deductive and rigorous microscopic ar¬ 
gument establishing its correctness is currently not avail¬ 
able [2], with the only exception of the integrable, non¬ 
interacting, case. Even then, the derivation is very non¬ 
trivial, and it has been provided only recently [13] using 
Discrete Holomorphicity (DH) methods. 

In this letter we present the first mathematical justi¬ 
fication of the quantum field theory description of non- 
integrable dimer models. We prove a theorem establish¬ 
ing the convergence, in the scaling limit, of the height 
function of model (1) to the massless Gaussian Free Field 
(GFF), in a suitable range of parameters. This is done by 
a new method, based on the combination of DH methods 
with Constructive Renormalization Group (CRG) tech¬ 
niques [14], which can be applied in a much wider con¬ 
text, including interacting dimers on different lattices and 
non-integrable deformations of Ising models. 

Given a dimer covering M, two faces of A centered at 
X and y and a path C'x-^.y from x to y with trivial wind¬ 
ing around the torus A, we define the height difference 
between x and y as 

hx — hy = (]lt,(M) ~ (2) 

b^Cx—¥y 
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where = + 1 / — 1 depending on whether Cx-^y crosses 
b with the white site on the right/left. Moreover, ]lb(M) 
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FIG. 1. A dimer configuration for L = 4 and the associated 
height function. The height of the central plaquette is con¬ 
ventionally set to 0. 

is equal to 1 if 5 is occupied by a dimer in M, and 0 
otherwise. A crucial property of the height function is 
that hx — hy is independent of the choice of Cx-s-y. The 
dimer eorrelation is given by (Ih; If,/), where is 

the truncated expectation with the weight in ( 1 ), and the 
two point height correlation is 

{ ihx ~ hy ) )= ^ ^ ' 0’6i 0'62 i ^62)- ( 3 ) 

{-y i-y 

Our main result is the following. 

Theorem. For A yf 0 suffieiently small, L —>■ 00, and 
m —> 0 , the height eorrelation for x 7 ^ y verifies: 

((hx-/iy)^) = ^^^log|x-y|-hi?(x-y), (4) 

with K{X) an analytie function such that K{0) = 1, 
and i?(x) a bounded remainder. The higher order trun¬ 
cated correlations are bounded uniformly in |x — y|. At 
large distances, the coarse graining ofhy^ converges to the 
Gaussian Free Field (GFF), in the sense that, if a € 
and f is a smooth, compactly supported function on 
with Jjg 2 f{u)du = 0 , one has 

/ f(u)f(v) log \u-v\dudv^ 

where s~^ represents the coarse-grain scale, to be sent to 
infinity after the thermodynamic limit. 

The choice of the specific interaction in (1) is just for 
illustrative purposes: the same result remains valid for 
generic finite range interactions, translationally and ro- 
tationally invariant. 

Eq.(5) can be re-read in a more evocative form: if xo is 
a smooth, compactly supported, probability distribution 
centered at the origin, X{(u) = Xo{u — f) is its translate, 
and ^xX{(ex), then choosing f = Xi - Xv 

in (5) we find 

^ {const.)\f - 


asymptotically as \^ — p] —)• 00. The left side is the 
coarse-grained “electric correlator”: our theorem proves 
its anomalous power law decay at large distances. 

An important step in proving the above result is the 
computation of the asymptotic behavior of the dimer cor¬ 
relation by CRG methods [14, 15]. In the limit L —> 00, 
if m yf 0 , it decays exponentially at large distances with 
rate (gaseous phase), with r]m{X) an ana¬ 

lytic function such that rim{0) =0. If m —> 0, it decays as 
a power law (liquid phase): e.g., if b, b' are both horizon¬ 
tal with b' — b = (xi,X 2 ) and z = Xi-\- ix 2 , then it decays 
polynomially, with critical exponent min{2, 2 -|-r 7 (A)}, and 
r]{X) = —(32/7r)A -|- O(A^) an analytic function of A, 

(1.; (6) 

Here K{X) is the same as in (4), K{X) is another an¬ 
alytic function such that K(0) = 1 , and h.o. indicates 
faster decaying terms at infinity. The above formula re¬ 
duces as A —> 0 to the one known by Kasteleyn’s exact 
solution. The main effect of the interaction is to produce 
an anomalous exponent in the second term, in agreement 
with the numerical simulations of [7]. Remarkably, there 
are no radiative corrections to the exponent of the first 
term. The model belongs to the same universality class 
as the XXZ chain, vertex models and Luttinger liquids. 

While the dimer characteristic function is a local ob¬ 
servable, the height differences are non-local “string” ob¬ 
servables, as apparent from (3). Even at A = 0, the com¬ 
putation of the height correlation is very subtle. Indeed, 
by inserting the A = 0 version of ( 6 ) into (3), one gets 
an apparently very singular expression: take e.g. x and 
y on the same horizontal line. In the large separation 
limit, the object of interest is formally proportional to 
// (m-vF 1 where ^,77 are the (suitably rescaled) horizon¬ 
tal coordinates of x and y. Such an integral requires a 
proper interpretation, because of its singularity at u = v, 
and the result depends on the specific ultraviolet reg¬ 
ularization. Of course, an “ad hoc” regularization can 
be chosen [16] in order to reproduce the expected result, 
but the problem remains of a general derivation, which 
can unambiguously return the correct exponents without 
any external bias. The problem was finally solved in [13], 
and the (I/tt^) factor in front of the logarithm in (4) at 
A = 0 was rigorously computed, by taking advantage of 
DH (lattice) methods. In the interacting case, the prob¬ 
lem is much more puzzling. In fact, in addition to the 
problem of the ultraviolet divergences affecting the com¬ 
putation of the (I/tt^) prefactor, the anomalous decay in 
( 6 ), once inserted into (3), may change the logarithmic 
growth into an anomalous growth. Our theorem proves 
that this is not the case: logarithmic fluctuations are ro¬ 
bust, stability being guaranteed by sophisticated cancel¬ 
lations arising from emerging chiral symmetry. Spurious 
ultraviolet divergences are avoided by using the irrelevant 
terms coming from the lattice: in this respect, the use of 
exact CRG methods (which, in contrast to field theoretic 
RG, takes the irrelevant terms into full account) is es- 
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sential. A detailed proof of our main theorem is rather 
technical and is given elsewhere [15]; below we explain 
its main ideas. 

Sketch of the proof. The first step consists in an ex¬ 
act rewriting of the finite volume/finite lattice gener¬ 
ating function of dimer correlations, 2(A), (defined so 
that (Ifci;--- ^tbk) = logZ(A))U^o, h label¬ 

ing the nearest neighbor bonds) as a finite Grassmann 
integral [15, Section 2]: 

Z(A) = ly2C0,r f ( 7 ) 

^ a.. Je,T 


Here are Grassmann variables, V is sum of monomials 
in A of order 4 or higher, B('ip,A) is a source term, 
sum of monomials in if and in A, 6/t G {0,1} label 
the boundary conditions for the Grassmann variables in 
the horizontal/vertical directions ( 0/1 corresponding to 
periodic/antiperiodic conditions), and Co,o = while 
= +1 otherwise. By cluster expansion methods, 
we prove that V and B are analytic in A. Pe^T-idtp) is a 
gaussian Grassmann integration with propagator g(x, y) 


1 _ji^(x-y) *sinfci-I-sinfc 2 + w(—1)^1 cosfci 

® 2 iA(k, to) 


( 8 ) 


where I?(k,TO) = to^ -I- (1 — TO^)(sinfci)^ -I- (sinA: 2 )^, and 
fci, ^2 are in (27r/L)Z or (27r/L)(Z -|- 1/2), depending on 
boundary conditions. 

If A = 0, then H = 0, in which case the integral is 
gaussian and can be computed exactly. When A 7 ^ 0 the 
integral is not gaussian, and it can be evaluated by a 
multiscale analysis using GRG methods [14]. We are in¬ 
terested in the case of to small or vanishing. As L —)■ 00 
and TO —>■ 0, the propagator in ( 8 ) becomes singular in 
correspondence of four momenta, namely pi = ( 0 , 0 ), 
P 2 = ( 7 r, 0 ), P 3 = ( 7 r, 7 r), p 4 = ( 0 , 7 r). Therefore, g(x,y) 
can be naturally written as the superposition of four 
terms, each of which is concentrated in momentum space 
around one of the singularities. Gorrespondingly, we de¬ 
compose the Grassmann held as: 


where V’x ,7 are Grassmann variables, often referred to 
as Majorana variables, since their effective action is a 
lattice regularization of the standard 2D Majorana ac¬ 
tion. Their propagator is block-diagonal, the helds with 
7 = 1,2 being independent of 7 = 3,4; the propagator 
G(x — y) of the 7 = 1,2 helds is the same as that of the 
7 = 3,4 helds, and reads (using the symbol / dk/(27r)^ 
as a shorthand for the discrete sum in ( 8 )): G(x) = 


1 

f dk 

X(k)e-**"^ 

f i sin ki + sin k 2 

im cos ki \ 

Z j 

(27r)2 

2D{k, to) 

—im cos fci 

i sin fci — sin k 2 J 


where x(k) is a smoothed characteristic function of the 
set max{|A:i|, |A: 2 |} < 7 r/ 2 , and Z = 1. To evaluate the 


Grassmann integral (7) we use (9) and write the propa¬ 
gator G(x) as sum of propagators living on momentum 
scales 2^,h < 0. After integrating the scales 0,..., h -I- 1, 
the (6*,r) contribution to Z{A) is rewritten as 

AAA) f 

Jd,T 

where Pzu,rah has propagator G*^^^(x), dehned in the 
same way as G(x), with Z replaced by Z^, to by to^ 
and x(k) by x?i(h), a (smoothed) characteristic function 
of the set |k| < (7r/2)2^. The effective potential is: 

yW(V>) = A^ E V'x,l'0x,2V’x,3V'x.4 Pir., 

X 

where ir. indicates the irrelevant terms (non-local quar- 
tic terms, and terms of order 6 or higher in ipY Remark¬ 
ably, the kernels of the irrelevant terms in are ana¬ 
lytic in X provided that jZ/i+i/Z/j —1|, |A/i| are sufficiently 
small, as long as \mh\ < 2^: the proof of this fact uses 
fermionic cluster expansion methods, including the use of 
Gram-Hadamard determinant bounds. Similarly, under 
the same assumptions, the effective source B^’^l is ana¬ 
lytic in A. Its structure is expressed most easily by using 
Dirac rather than Majorana helds: the former are dehned 

^x,l ■“ T Ax,3)1 ^x, —1 ■“ ^72 ^ Ax,4), 

and they are referred to as Dirac variables, because their 
action is the lattice analogue of that of 2D Dirac helds. 
In terms the effective source reads: 

yip y(2) 

bW(A = + -^F2(rP,3) + ir. , 

where ir. are the irrelevant terms (non local, or of 
higher order in A or as compared to Aipip)- More¬ 
over, denoting Jx,i = <7(x,x-i-ei) with Jh = e^'’ — 1: 
Fi = 2Ex. <^=±(-l)''('7x.i + and F 2 = 

2Ex.^c;=± Sum¬ 

marizing, the effective theory on scale h has the same 
structure as a theory of interacting 2D lattice Dirac 
fermions with a wave function renormalization Zh, an 
effective mass irih, an effective coupling \h, and effective 
source couplings zl^^\ zj^^K It is completely analogous 
to that obtained in the multiscale analysis of the 8 Ver¬ 
tex, Ashkin-Teller, XXZ, or Luttinger liquid models [14]: 
the only differences have to be found in the oscillating 
factors appearing in the dehnition of Pi, F 2 and in the 
specihc structure of the irrelevant terms. The how equa¬ 
tion for the effective couplings of all these models is the 
same, up to irrelevant contributions, which are exponen¬ 
tially negligible in the infrared limit. Therefore, Xh ap¬ 
proaches exponentially, as /i —)■ — 00 , a line of hxed points: 
A_oo(A) = -32A(1 -b 0(A)). Moreover, Zh ~ 

Z^^ ~ ^ ,^ 2 Um(>dh^ where ~ means that 

the ratio of the two sides is bounded from above and be¬ 
low by two universal positive constants, uniformly in h. 
Remarkably, using the emergent chiral gauge symmetry 
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of the theory, we find that 77 = r/i, which implies the 
robustness (exact non-renormalization) of the exponent 
2 in the first term of ( 6 ). The integration goes on until 
nih ~ 2^, at which point the Dirac field is massive and 
can be integrated in one step. If m —>■ 0 and L —> 00 , the 
integration has no infrared cutoff. 

In order to evaluate the height fluctuations, we use 
the path-independence of the height difference, which 
is a (weak) instance of DH. We proceed as in [13](c). 
Consider e.g. the height variance: in the right side 
of (3) we deform the two paths along which bi and 
62 are summed over, in such a way that they are “as 
much separated as possible”, as in Fig.2. In the vicin- 



FIG. 2. A schematic view of the paths along which bi, 62 are 
summed over, to be called CxXy and C^Xy 


ity of X and y, the two paths are lattice approxima¬ 
tions of straight lines, departing from and arriving at 
the points x, y in different directions. After the path 
deformation, we replace the dimer correlation in the 
right side of (3) by its asymptotic expression ( 6 ) (and 
its analogues in the cases that 6 , b' have different ori¬ 
entations). The h.o. terms contribute a finite constant, 
uniformly in jx — yj. The contribution to (3) from the 
term with decay exponent 2 (let us call it Conff) reads: 


Conto = — 


_ K(X) 


27r2 


E 


bi GC: 


(1), Efc^eca. 


Zbi are the representatives in complex coordinates of the 
centers of the bonds bi, and are the oriented elemen- 
tary path elements of CxAy crossing bi, expressed in com¬ 
plex coordinates. Note that no oscillatory factor appears 
in Cont 2 ' the factors abi<Jb 2 in (3) compensate exactly 
the oscillatory factor of the term under consideration in 
the dimer correlation. Cont 2 is the Riemann approxi¬ 
mation to -'|:^Re dz dw , where 71 and 72 

are two completely disjoint complex paths (this is what 
makes the integral non-singular!) going from Zx = - 2 ^( 1 ) 

to Zy = z^(i), and from = z^( 2 ) to Zy = z^( 2 ), where bX^ 

(i) ^ (i) 

and by are the first and last bonds of CiAy. Its value 

is ^Re log \y_X\{XXX] j which is the same as (4) up 
to a bounded error. Finally, consider the contribution 
to (3) from the term with exponent 2 -|- 77 : in this case 
the factors do not compensate exactly with the 

oscillatory signs in the dimer correlation; the left-over 
oscillations act, after summation along the paths, as dis¬ 
crete derivative, which effectively makes this term decay 
faster, thus making its contribution to (3) finite, uni¬ 
formly in jx — yj. Similar considerations apply to higher 
order cumulants, and (5) follows as a corollary. 

In conclusion, we presented a rigorous microscopic 
derivation of massless gaussian free field behavior of the 
height field of a non integrable interacting dimer model. 
Our method combines constructive field theory tech¬ 
niques with discrete holomorphicity ideas, which are used 
for the first time in a unified way to analyze a non-local 
fermionic observable. The method can be applied to sev¬ 
eral other non-integrable 2D critical theories and we ex¬ 
pect it to be capable, in perspective, of rigorously proving 
conformal invariance of the scaling limit. 
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